Interpolative Approach for Solving Quantum Impurity Model Based on the 
Slave— Boson Mean— Field Approximation 
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A rational representation for the self-energy is explored to interpolate the solution of the Anderson 
impurity model in general orbitally degenerate case. Several constrains such as the Friedel's sum rule, 
high-frequency moments and the value of quasiparticle residue are used to establish the equations 
for the coefficients of the interpolation. We test two fast techniques, the slave-boson mean-field 
and the Hubbard I approximation to determine the coefficients. The obtained self-energies are 
compared with the results of numerically exact Quantum Monte Carlo method. We find that using 
the slave-boson mean-field approach we can construct an accurate self-energy for all frequencies 
via the proposed interpolation procedure. 

PACS numbers: 71.10.-w, 71.27. +a, 71.30.+h 



I. INTRODUCTION 

There has been recent progress in understanding 
physics of strongly correlated electronic systems and their 
electronic structure near a localization-delocalization 
transition through the development of dynamical mean- 
field theory (DMFT) 1 . Merging this computation- 
ally tractable many-body technique with realistic local- 
density-approximation (LDA)2 based electronic struc- 
ture calculations of strongly correlated solids is promis- 
ing due to its simplicity and correctness in both band 
and atomic limits. At present, much effort is being 
made in this direction including the developments of a 
LDA+DMFT method 3 -, LDA++ approach^, combined 
GW and DMFT theorji as well as applications to vari- 
ous systems such as Lai-xSr^TiOg^, V 2 0^, Fe and Ni&, 
Ce^, PuiSiii, and many others. For a review, see Ref. 

03 

The development of ah initio DMFT scheme requires 
fast methods and algorithms to solve the Anderson im- 
purity model! 3 in general multiorbital case. Present 
techniques based on either non-crossing approximation 
(NCA) or iterative perturbation theory (IPT) are unable 
to provide the solution to that problem due to a limited 
number of regimes where these methods can be applied 1 . 
The Quantum Monte Carlo (QMC) technique-*!* is very 
accurate and can cope with multiorbital situation but 
not with multiplet interactions. Also its applicability 
so far has been limited either to a small number of or- 
bitals or to unphysically large temperatures due to its 
computational cost. Recently some progress has been 
achieved using impurity solvers that improve upon the 
NCA approximationiSiiStil, but it has not been possi- 
ble to retrieve Fermi liquid behavior at very low tem- 
peratures with these methods in the orbitally degenerate 
case. 

In this paper we explore the possibility to use interpo- 



lation for the self energy of the quantum impurity model 
in general multiorbital situation. We do not attempt to 
develop an alternative method for solving the impurity 
problem, but follow the ideology of LDA theory where 
approximations were designed by analytical fits"* to the 
Quantum Monte Carlo simulations for homogeneous elec- 
tron gasi^. Numerically very expensive QMC calcula- 
tions for the impurity model display smooth self-energies 
at imaginary frequencies for a wide range of interactions 
and dopings, and it is therefore tempting to design such 
an interpolation. We also keep in mind that for many 
applications a high precision in reproducing the self- 
energies may not be required. One of such applications 
is, for example, the calculation of the total energja^ifliii 
which, as well known from LDA based experience, may 
not be so sensitive to the details of the one-electron spec- 
tra. As a result, we expect that even crude evaluations 
of the self-energy shapes on imaginary frequency axis 
may be sufficient for solving many realistic total energy 
problems. Another point is a computational efficiency. 
Bringing full self-consistent loops with respect to charge 
densities"* and other spectral functions require many it- 
erations towards the convergency which may not need 
too accurate frequency resolutions at every step. This 
view should, of course, be contrasted with calculation 
of properties such as the low-energy spectroscopy and 
especially transport where delicate distribution of spec- 
tral weight at low energy, namely the imaginary part of 
the analytically continued self-energy, needs to be com- 
puted with much greater precision. Here, extensions of 
the interpolative algorithms should be implemented and 
its beyond the scope of the present work. 

We can achieve a fast interpolative algortihm for the 
self-energy utilizing a rational representation. The coef- 
ficients in this interpolation can be found by forcing the 
self-energy to obey several limits and constrains. For ex- 
ample, if infinite frequency (Hartree-Fock) limit, high- 
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frequency moments of the self-energy, low-frequency 
mass renormalization, number of particles as well as the 
value of the self-energy at zero frequency are known from 
independent calculation, the set of interpolating coeffi- 
cients is well defined. In this work, we explore the slave- 
boson Gutzwiller approacbS&SLS^iSii and the Hubbard I 
approximation^ to determine these coefficients. Com- 
paring interpolating results with more accurate calcula- 
tions using the Quantum Monte Carlo method we find 
that the slave-boson approach predicts the parameters 
of interpolation with a good accuracy while the Hubbard 
I method fails in a number of regimes. Thus, interpo- 
lation for the self-energy at all frequencies can be ob- 
tained from the method designed to work exclusively at 
low-frequency limit . 

Our paper is organized as follows. In Section II we dis- 
cuss rational interpolation for the self-energy and list the 
constraints. In Section III we discuss methods for solv- 
ing Anderson impurity model based on the slave-boson 
Gutzwiller and the Hubbard I approximations which can 
be used to find these constraints. We also present our 
numerical comparisons with the QMC data between var- 
ious quantities such as the number of particles and the 
quasiparticle residue. In Section IV we discuss the inter- 
polational approach and compare the self-energies and 
the Green functions with the QMC results. Section V is 
the conclusion. 



II. INTERPOL ATIVE APPROACH 

To be specific, we concentrate on the Anderson impu- 
rity Hamiltonian 

JV ^ N ^ 
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describing the interaction of the impurity levels e/« 
with bands of conduction electrons E ka via hybridization 
V^(k). U a p is the Coulomb repulsion between different 
orbitals in the /-band. Inspired by the success of the 
iterative perturbation theory!, in order to solve the An- 
derson impurity model in general multiorbital case, we 
use a rational intcrpolative formula for the self energy 
This can be either encoded into a simple form of pole 
expansion 



T, a (iu) = E Q (ioo) + 2J 



iuj - PS 



(2) 



or, alternatively, in a form of continuous fraction 
expansion 



Yi a {iw) = E a (ioo) + 



iuj — B a 



(3) 



The coefficients W", P™, in Eq. © or the coefficients 
A a , B a ,C a , ... in Eq. © are to be determined to- 
gether with the Hartree Fock value of the self-energy 
E a (ioo). In principle, both representations are equiva- 
lent and there is a well-defined non-linear relationship 
between the parameters entering J2J) and J2J. 

Our basic assumption is that only a few poles in the 
rational representation J2J) or a few coefficients in con- 
tinuous fraction expansion J3J is necessary to reproduce 
an overall frequency dependence of the self-energy nec- 
essary for the accuracies of the total energy calculations. 
Extensive experience gained from solving Hubbard and 
periodic Anderson model within DMFT at various ratios 
of the on-site Coulomb interaction U to the bandwidth 
W shows the appearance of lower and upper Hubbard 
bands as well as renormalized quasiparticle peak in the 
spectrum of one-electron excitations^. To describe such 
three-peak structure we expect that either two- or at 
most three-pole formulae for the self-energy will suffice. 
Thus, four (six) unknown parameters need to be deter- 
mined when using two (three) pole interpolating formu- 
lae. 

In order to fix the coefficients we can explore several 
constrains for the self-energy: 

a) Hartree-Fock value. In the limit iuj — > ioo the self- 
energy takes its Hartree-Fock form 



13 



(4) 



b) Moments. The subleading high frequency behavior 
of the self energy can be used to determine some addi- 
tional coefficients. In the case of the three pole approx- 
imation, we determine coefficients A a ,B a and C a from 
the first three moments of the self-energy Y$ , , Sa' 
as follows 



A - E (1) 

"-a — ) 



B a = 



v( 3 )vW (V^>\ 



(2h2 



(•5) 



(0) 



(7) 



■J-D a 



while in the two pole approximation only coefficient A a 
should be determined from the high frequency moment, 
as stated in Eq. (JjjJ above. The self-energy moments can 
be expressed by the occupancies and correlation func- 
tions that can be calculated within an approximation like 
slave-boson Gutzwiller or Hubard I. The expressions for 
the self-energy moments are worked out explicitely in the 
Appendix A. 
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c) Mass Renormalization. If independent determina- 
tion of the quasiparticle residue, Z a , exists, the following 
relationship holds 



-g- U=o- 1 - Z a . 



(8) 



d) Number of particles. Sum over all Matsubara fre- 
quencies of the Green's function gives number of parti- 
cles, i.e. 



where 



G fa (iuj) 



- Cf a - A a (iu) - T, a (iu>) 



(9) 



(10) 



defines the impurity Green function and A a (iuj n ) is the 
hybridization function. 
e) Friedel sum rule. 

This is a relation between the total density and the real 
part of the self-energy at zero frequency. It determines 
zero-frequency value of the self-energy 



1 1 fe fa + 3?£ a (i0+) + 3fJA Q (i0+) 

na = 2 + ^ arctg { 



-\-ioc 



dA a {z) ^ z0+ 
dz 



(11) 



Formally, some of the constrains [(c) and (e)] hold for 
zero temperature only but we expect no significant devi- 
ations as long as we stay at low temperatures. 

We thus see that the interpolational scheme is defined 
completely once a prescription for obtaining parameters 
such as Z Q , n a and Sq™^ is given. For this purpose we 
will test two popular methods: slave-boson Gutzwillcr 
method** as described by Kotliar and Ruckenstein** and 
the well-known Hubbard I approximation**. We com- 
pare these results against more accurate but computa- 
tionally demanding Quantum Monte Carle 1 calculations. 

Notice that once the parameters such as Z a are com- 
puted from a given approximate method, some of the 
quantities such as the total number of particles, n a , and 
the value of the self-energy at zero frequency, £ Q (z0), 
can be computed fully self-consistently. They can be 
compared with their non-self-consistent values. If the 
approximate scheme already provides a good approxi- 
mation for n a and satisfies the Friedel sum rule, the 
self-consistency can be avoided hence accelerating the 
calculation. Indeed we found that inclusion of the self- 
consistency improves the Gutzwiller results in the vicin- 
ity of the Mott transition only where the Hubbard bands 
get increasingly important. 

We now give the description of our approximate algo- 
rithms and then present the comparison with the QMC 
calculations. 



III. METHODS FOR SOLVING IMPURITY 
MODEL 

A. Quantum Monte Carlo Method 

The Quantum Monte Carlo method is a powerful and 
manifestly not perturbative approach in either interac- 
tion U or the bandwidth W. In the QMC method one in- 
troduces a Hubbard-Stratonovich field and averages over 
it using the Monte Carlo sampling. This is a controlled 
approximation using different expansion parameter, the 
size of the mesh for the imaginary time discretization. 
Unfortunately it is computationally very expensive as 
the number of time slices and the number of Hubbard- 
Stratonovich fields increases. Extensive description of 
this method can be found in Ref. 01 We will use this 
method to benchmark our calculations with approximate 
but much faster algorithms described below. 



B. Slave-Boson Approach 

A fast approach to solve a general impurity problem 
is the slave-boson method**^***. At the mean field 
level, it gives the results similar to the famous Gutzwiller 
approximation**. However, it is improvable by perform- 
ing fluctuations around the saddle point. This approach 
is accurate as it has been shown recently to give the ex- 
act critical value of U in the large degeneracy limit at 
half-filling**. 

The main idea is to rewrite atomic states consisting 
of n electrons |7i,...,7„), < n < N with help of a 
set of slave-bosons {ipZ 1 '' " ,7 "}. In the following, we as- 
sume SU(N) symmetric case, i.e., equivalence between 
different states |7i,...,7n) for fixed n. Formulae corre- 
sponding to a more general crystal-field case are given in 
Appendix B. The creation operator of a physical electron 
is expressed via slave particles in the standard manner**. 
In order to recover the correct non-interacting limit at 
the mean-field level, the Bose fields ip n can be considered 
as classical values found from minimizing a Lagrangian 
L{ip n } corresponding to the Hamiltonian Two La- 
grange multipliers A and A should be introduced in this 
way, which correspond to the following two constrains: 



N 



N 



(12) 



3 M0+ = n. (13) 



The physical meaning of the first constrain is that the 
sum of probabilities to find atom in any state is equal to 
one, and the second constrain gives the mean number of 
electrons coinciding with that found from G g (iuj) 



A combinatorial factor = 



Nl 



A - b 2 A(lu)Y . ~ ^ux^uux^x x™ u„ - ra!(j v_ ra) , 

arrives due to assumed equivalence of all states with n 
electrons. 
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Minimization of L{ip n } with respect to ip n leads us to 
the following set of equations to determine the quantities 

ipn- 

[E n +K~n\]^ n +nbT^k{iu)G g {iu)[LR4i n _ l +4i n bL 2 ] + 



+{N - n)bT A(iu)G g (iuj)[R 2 bip n + LRip n+1 ] = 0, 

(14) 

where b — ^ E n =i C^-iVn'/'n-i, determines the 
mass renormalization, and the coefficients L = (1 — 

ELiC^lO- 1 ' 8 . * = (1 - ElaC^r 1 / 2 are 
normalization constants as in Refs. 121123 . -E„ = e/n + 
Un(n — l)/2 is the total energy of the atom with n elec- 
trons within SU (TV) approximation. 

Eq. 1(14(1 . along with the constrains If 2J1 . (| 1 3(1 consti- 
tute a set of non-linear equations which have to be solved 
iteratively. In practice, we consider Eq. (f 1 4|> as an eigen- 
value problem with A being the eigenvalue and ipn being 
the eigenvectors of the matrix. The physical root corre- 
sponds to the lowest eigenvalue of A which gives a set of 
ip n determining the mass renormalization Z = b 2 . Since 
the matrix to be diagonalized depends non-linearly on ip n 
via the parameters L,R, and b and also on A, the solu- 
tion of the whole problem assumes the self-consistency: 
(i) we build an initial approximation to ip n (for exam- 
ple the Hartree-Fock solution) and fix some A, (ii) we 
solve eigenvalue problem and find new normalized ip n , 
(iii) we mix new ip n with the old ones using the Broyden 
method^ and build new L, R, and b. Steps (ii) and (iii) 
are repeated until the self-consistency with respect to ip n 
is reached. During the iterations we also vary A to obey 
the constrains. The described procedure provides a sta- 
ble computational algorithm for solving AIM and gives 
us an access to the low-frequency Green's function and 
the self-energy of the problem via knowledge of the slope 
of QT,(iu)) and the value 3fJ£(0) at zero frequency. 

The described slave-boson method gives the following 
expression for the self-energy: 

= (l-b- 2 )iu}-e f + \b- 2 . (15) 

The impurity Green function Gf(iuj) in this limit is given 
by the expression 

G f {iuj) = b 2 G g {iuj). (16) 

As an illustration, we now give the solution of Eq. 1(14(1 
for non-degenerate case (N = 2) and at the particle-hole 
symmetry point with ei mp — fi = —^(N — I). Con- 
sider a dynamical mean field theory for the Hubbard 
model which reduces the problem to solving the impurity 
model subject to the self-consistency condition with re- 
spect to A(iuj). Starting with the semicircular density of 
states (DOS), the self-consistency condition is given by 
A(iw) = t 2 Gf(iuj) with parameter t being the quarter 
of bandwidth. We obtain the following simplifications: 



L = R = y/2, A = 0, Vo = ip2, b = 4^iV>2 and G g (iuj) = 
[lev - ft^jH]" 1 . The sum TJ2 iu A{ioj)G g {ioj) ap- 
peared in Eq. I(f 4|) scales as 2ta with the constant a 
being the characteristic of a particular density of states 
and approximately equal to -0.2 in the semicircular DOS 
case. The self-consistent solution of Eq. ((f4(l is there- 
fore possible and simply gives ip 2 = 12 g tQ , + \ ■ The Mott 
transition occurs when no sites with double occupancies 
can be found, i.e. when ipo = ip2 = 0. A critical value of 
U c = 32t\a\ = 8W\a\ with W being the bandwidth. Fol- 
ey ps —0.2, this gives U c « 1.6W and reproduces the exact 
result for U c — 1.49VF within a few percent accuracy. As 
degeneracy increases, critical U is shifted towards higher 
values. From numerical calculations we obtained the fol- 
lowing values of the critical interactions in the half-filled 
case U c « 3W for N = 6 (p-level), U c w 4.5W for N = 10 
(d-level), and U c w 6W for N = 1 4 (/-level). 

Density-density correlation function (n a np) (where 
a / |3) for local states with n electrons is proportional 
to the number of pairs formed by n particles CQ/C^. 
Since the probability for n electrons to be occupied is 
given by: P n = ip 2 Cn, the physical density-density cor- 
relator can be deduced from: (nn) = ^ZnC^jC^Pn- 
Similarly, the triple occupancy can be calculated from: 
(nnn)=J2 n C2/Ci r P n . 

Let us now check the accuracy of this method by com- 
paring it's results with the QMC data. We consider the 
two-band Hubbard model within St/ (4) orbitally degen- 
erate case. Hybridization A(ioj) — (k) / (iuj — Sk) 
satisfies the DMFT self-consistency condition of the 
Hubbard model on a Bethe lattice 

A(iu) =t 2 G(iuj), (17) 

where t is a quarter of the bandwidth W . The Coulomb 
interaction is chosen to be U = Ad (d is the half- 
bandwidth) which is sufficiently large to open the Mott 
gap at integer fillings. All calculations are done for the 
temperature T = 1 /16rf. 

We first compare the average number of electrons 
vs. chemical potential determined from the slave bosons 
which is plotted in Fig. ^a). This quantity is sensitive 
to the low-frequency part of the Green function which 
should be described well by the present method. We see 
that it reproduces the QMC data with a very high accu- 
racy. The quasiparticle residue Z = 1/(1— ^gffi" )lto;->o 
versus filling n is plotted in Fig. db). The slave-boson 
method gives the Fermi liquid and provides a good es- 
timate for the quasiparticle residue. Indeed, the agree- 
ment with the QMC results is satisfactory. Also, recent 
findings 25 suggest that at half-filling Z deduced from 
slave bosons becomes exact when N — > oo. Still the dis- 
crepancy of the order of 30% can be found. In fact, we 
have performed several additional calculations for other 
degeneracies (N = 2 and 6) and for various parameters 
regimes. The trend to underestimate mass renormaliza- 
tion by about 30% can be seen in almost all cases (it 
disappears only when U approaches zero). This empiri- 
cal finding is useful as we can use the reduced values of 
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n 

FIG. 1: Comparison between the slave boson mean field and 
the QMC calculation for (a) concentration versus chemical 
potential, (b) dependence of the spectral weight Z on concen- 
tration, and (c) density-density correlation function, {n a n a i) 
versus filling, n, in the two-band Hubbard model. 

Z needed for our interpolating formula. 

Fig. ^c) shows the density-density correlation func- 
tion (n a n a i) as a function of average occupation n. The 
discrepancy is most pronounced for fillings h < 1 [sec the 
inset of Fig. ^c)] where the absolute values of (n a n a i) 
are rather small. Although our slave-boson technique 
captures only the quasiparticle peak, it gives the correla- 
tion function in reasonable agreement with the QMC for 
dopings not too close to the Mott transition. 

C. Hubbard I Approximation 

Now we turn to the Hubbard I approximations* which 
is closely related to the moments expansion method 27 . 
One introduces Hubbard operators and the correspond- 
ing Green function G nm (iu) defined for these Hubbard 
operators. Within SU(N) we are summing over certain 
set of Hubbard operators which induce transitions be- 
tween n-times and n + 1-times occupied local level. The 
impurity Green function is 

Gf(iu) =J2 G nm(iu), (18) 

nm 

where the matrix [G nm (iw)] _1 = [G°*„(io;)] _1 — 
A(iu) 7 n,m — 0,N — 1. Gf* m (iu) is the atomic Green 
function 

IU + (J, - E n+ i + E n 



The coefficients X n are the probabilities to find atom 
with n electrons similar to the coefficients i/> 2 introduced 
above. They can be also used to find the averages (ran) 
in the same way. These numbers are normalized to 
unity, Et Cn*» = T^C^- 1 (X n +X n+1 ) = 1, 
and are expressed via diagonal elements of G nm (iu) 

as follows: X n = -T^ G^He"^/^" 1 . The 
mean number of electrons can be measured as follows: 
n = En=o nC^X n or as: h — TN £ iw G f (iu)e^ 0+ . 

Note that when A(iu) = 0, Gf(iu) is reduced to 
J2 nm 4Hi i- e - Hubbard I reproduces the atomic 
limit. Setting U = gives Gf(iu>) = [iuj + /i — ef — 
A(iw)] -1 , which is the correct band limit. Unfortu- 
nately, at half-filling this limit has a pathology con- 
nected to the instability towards Mott transition at any 
interaction strength U. To see this, we consider a dy- 
namical mean-field theory for the Hubbard model. Us- 
ing semicircular density of states, we obtain Gf{iuj) — 
[1 — t 2 Gf (iu) G at (iu)]~ 1 G a t(iu) and conclude that for 
any small U the system opens a pathological gap in the 
spectrum. Clearly, using Hubbard I only, the behavior 
of the Green function at iu — > cannot be reproduced. 
This emphasizes the importance of using the slave-boson 
treatment at small frequencies. 

Let us now check the accuracy of this method against 
QMC. We again consider the two-band Hubbard model 
within 5 f C/(4) symmetry as in the case of the slave- 
boson method described earlier. The n{jj,) is plotted in 
Fig. [2Ja) . The Hubbard I approximation does not give 
satisfactory agreement with the QMC data because it 
misses the correct behavior at low frequencies. The com- 
parisons for Z(n) is plotted in Fig. Hfb). These num- 
bers are extracted from the low-frequency slope of the 
imaginary part of the self-energy on the imaginary axis. 
Surprizingly, the Hubbard I Z's show a relatively good 
behavior when comparing to the QMC method. How- 
ever, the patology of this approximation at half-filling 
would predict Z — for any U, which is a serious warn- 
ing not to use it for extracting the quasiparticle weight. 
Fig. |5fc) shows (nn) as a function of average occupa- 
tion n. As this quantity is directly related to the high- 
frequency expansion one may expect a better accuracy 
here. However, comparing Fig. j^c) and Fig. \Rc), it 
is clear that the slave boson method gives more accu- 
rate double occupancy. This is due to the fact that the 
density matrix obtained by the slave boson method is of 
higher quality than the one obtained from the Hubbard 
I approximation. 

The results of these calculations suggest that the slave- 
boson Gutzwiller method shows a good accuracy in com- 
parison with the QMC. Therefore it can be used to de- 
termine the unknown coefficients in the intcrpolational 
form of the self-energy, Eqs. @ or J2J). At the same 
time the use of the Hubbard approximation should be 
avoided. Interestingly, while more sophisticated QMC 
approach captures both the quasiparticle peak and the 
Hubbard bands this is not the case for the slave-boson 
mean field method. To obtain the Hubbard bands in this 




ttQMC 
■° Hubbard 



n 



FIG. 2: Comparison between the Hubbard I and the QMC 
calculation for (a) concentration versus chemical potential, 
(b) dependence of the spectral weight Z on concentration, 
and (c) density-density correlation function, (n a n a i) versus 
filling, n, in the two-band Hubbard model. 



method fluctuations need to be computed, which would 
be very tedious in a general multiorbital situation. How- 
ever the slave-boson method delivers many parameters 
in a good agreement with the QMC results, and hence it 
can be used to obtain the parameters that enter the ra- 
tional approximation. In this way slave-boson approach 
generates a better Green function with quasiparticles and 
Hubbard bands at no extra computational cost. 



IV. RESULTS OF THE INTERPOLATIVE 
SCHEME 

We now turn to the comparison of the self-energie ob- 
tained using the formula @ and the corresponding Green 
functions after the formula (|1U[) against the predictions 
of the Quantum Monte Carlo method. We have studied 
the results of the two-pole and the three-pole interpo- 
lations which require the determination of either 4 or 6 
coefficients. The interpolational formula with two-poles 
needs finding the values of S(0), Z and the first high- 
frequency moment of the self-energy E*- 1 ). The condition 
that the total number of electrons found from the Green 
function (|10l) is the same as predicted by the slave-boson 
method can be used to fix the fourth coefficient. To ex- 
plore the three-pole interpolation, two more conditions 
are required. We have used additional high frequency 
moments E^ 2 \ and T,^ which can be determined if one 
computes correlational functions (ran), (nnn), (nnnn). 
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FIG. 3: Comparison between real and imaginary parts of the 
Green function and the self-energy obtained from the inter- 
poaltive method and the Quantum Monte Carlo calculation 
for the two-band Hubbard model at filling n — 0.8. 

We perform the comparison of the interpolational 
scheme against QMC for two characteristic values of the 
concentrations: n — 0.8 representing the case n < 1 and 
n = 1.2 representing the case n > 1. Both cases are 
chosen to be close to the Mott transition. This is a hard 
test for our approach as we have found that the pre- 
dictions of the interpolative scheme and the QMC data 
are practically indistinguishable in the metallic region 
far from the Mott transition (U < W). The two-band 



0.6 



,0.4 



O 

a 



0.2 



1 


l 


i 

\ — QMC r 


i 






\ 2 pole p 




; | 




1 — 3 pole 1 








1 n f =0-8 ! 




J, 




1 . 1 


i ! 


4 -2 





2 


4 ( 



FIG. 4: Comparison between the one-electron densities of 
states obtained from interpolative formula and the Quantum 
Monte Carlo calculationfor the two-band Hubbard model at 
filling n = 0.8. 

Fig- OH shows the comparison of the real and imaginary 
parts of the self-energy and the Green function obtained 
by the two- and three-pole interpolations with the QMC 
results. Both interpolational schemes give pretty good 
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agreement at all frequencies. These results are obtained 
with the values of mass renormalizations Z reduced by 
30% which gives a better agreement with Zqmc as we 
discussed in Section III. The results with bare Z's de- 
duced from the slave-boson method are also found to 
be quite accurate but not as good as we advertise in 
Fig. |3 We have also performed the comparison for the 
one-electron density of states which are shown in Fig. 0] 
One can see the appearance of the lower and upper Hub- 
bard bands as well as renormalized quasiparticle peak. 
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FIG. 5: Comparison between real and imaginary parts of the 
Green function and the self-energy obtained from the inter- 
poaltive method and the Quantum Monte Carlo calculation 
for the two-band Hubbard model at filling n — 1.23. 

In Fig. \E\ we show our results for the case n — 1.2 
using the two- and three-pole interpolational formula 
©. As far as the two-pole scheme is concerned, the 
overall agreement is satisfactory but a discrepancy exists 
in the region of intermediate imaginary frequencies. It 
is directly related to the misplaced upper Hubard band 
as can be seen in Fig. [fj] This discrepancy disappears 
in the three-pole interpolation as is evident from the 
Fig. [5] Comparison of the one-electron density of states 
is shown in Fig. Clearly, the two-pole interpolation 
is inadequate in this regime, while the three-pole inter- 
polating scheme gives correct position of Hubbard bands 
and quasiparticle peak, all basic features of the Anderson 
impurity problem. 

Finally, we would like to comment on the renormaliza- 
tion of the quasiparticle residue Z beyond the slave boson 
mean-field result. While we found empirically that the 
30% renormalization significantly improves our fits to the 
QMC data, over a broad region of parameters, this can- 
not hold everywhere since Z should go to 1 when U ap- 
proaches zero. The reduction of Z beyond the mean-field 
results, in the regions where the correlations are strong, 
can be thought to arise from fluctuation corrections of 



the slave boson around the saddle point value. A com- 
putationally efficient way to evaluate these corrections, 
which will also give rise to the quasiparticle lifetime not 
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CO 

FIG. 6: Comparison between the one-electron densities of 
states obtained from interpolative formula and the Quantum 
Monte Carlo calculationfor the two-band Hubbard model at 
filling n = 1.23. 



V. CONCLUSION 

To summarize, this paper shows the possibility to ex- 
tract parameters for the interpolating self-energy of the 
Anderson impurity model such as the self-energy at zero 
frequency, the mass renormalization, the density and the 
high frequency moments using computationally very ef- 
ficient slave-boson mean field approach. These data are 
then used to reconstruct the self-energy at all frequen- 
cies using a rational interpolation formula. We find that 
three-pole interpolation predicts the self-energy quite ac- 
curately as compared to the results of more sophisticated 
approaches such the Quantum Monte Carlo method. 

As a general conclusion, we have thus found that inter- 
polative formulae with several poles are sufficient to fit 
the self-energies in a broad range of the parameters of the 
quantum impurity problem and the slave-boson mean- 
field method is capable to deliver all necessary ingredi- 
ents for this algorithm. The approach can be used for the 
evaluation of LDA+DMFT total energies and is a very 
direct way to obtain spectra with a three-peak structure, 
characteristic of the strongly correlated regime, with very 
little computational effort. 
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VI. APPENDIX A 



High-frequency moment expansions for the Green function and the self-energy can be worked out: 



G a (iu>) 



r W r>(2) r (3) 



(20) 



where the moments of the Green function are given by the following formulae 
G« = 1 



00' k 

G i 4) = <4a + 3e 2 a U a 0(ri0) + 3e /a U a 0U a 0> {n07ip,) + U a pU a p>U a p>i{npTip>npn 
p 00' pp'P" 

+ ]T |K(k)| 2 (£ k a + 2e /a + 2^2 U a p{np)) + ]T C/^(k)((2n a - 1) /Jq^). 

k /3 k,/3 



(21) 
(22) 

(23) 



(24) 



The expansion for the self-energy is obtained as follows 



E a (iu) = + /i — A a (iu>) 



' r W r (2) rt(3) 

W 2 



Er»(ioo) 



v (i) v (2) v (3) 



(25) 



where the self-energy moments are given by 



4 15 

S (2) 



Gi 3) - [Gi 2 )] 2 -El^! 2 = E UapUvpifonp,) - (np)(np,)), 
k 00' 

G$ - 2G^G^ + [G^] 3 - £ |V Q (k)| 2 £ ka = 

k 

^ U a 0U a 0,U a 0„ {{710710,710,,) - 2{npri0,){n0„) + {np) (np>) (npi>)) 
pp'P" 

e !c t ^U a pU a p l {{npnp l ) - {np){np>)) + ^C/^V^(k)((2n Q - 1)^^), 



g! 5 )-2g1 4 )g1 2 ) 



GL 3 f + 3G( 3 ) [Gi 2 f - [G^f-^ly^k)! 2 ^ 



r 2 

J ka 



2J ^a/Jl ^afe U a 3 U a 4 { {np 1 Tl0 2 Tl0 3 71 4 ) - 2 {ri0 1 71 p 2 71 p 3 ) {np l ) - {np^p^ {flp 3 Tl0 4 

0102 0304. 

+ 3{n0 1 n0 2 ){n0 3 ){n0 i ) - {np l ){np 2 ){np s ){np /k )) 
2e a 2J U a 1 U a 2 U a 3 {{n0 1 n0 2 n0 3 ) ~2{n0 1 7i0 2 ){n0 3 ) + {7i0 1 ){n0 2 ){n0 3 )) 



010203 



€ l Yl U a0i U <*02 i( n 0i n 02) - { n 0l){ n 02)) +A/^ 3) . 

0102 



(26) 



(27) 



(28) 



The nonlocal part of the self-energy moment, which explicitly depends on hybridization matrix V and vanishes in the 
atomic limit, is explicitly quoted only for the second moment (last term in Eq. |23). In the third order, the nonlocal 
term is denoted by Aa and has not yet been evaluated. In practical calculations, we neglected it's contribution. 
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Within SU(N) symetric case the formulae are simplified. First, introduce the definitions: 



N 

n 



<»> = E # ( M ) 



n=l 

N 



fin 

= E # ( 3 °) 



("»»> = (31) 
- = ^E C ™-2 2 ^-^« (32) 

n— 1 

where ^ is probability for n times occupancy (with n electrons on the local level). Then 

Y J U a p{n l3 ) = U{N-l){n), (33) 



E U a pU a p,{n p n fi ,) = U 2 [(N - 1)(N - 2){nn) + (N - l)<n>] , (34) 

00' 

E U af3 U a0 ,U af} „(n n m ,) = U 3 [(N - 1)(N - 2)(N - 3)(nnn) + 3(N - 1)(N - 2){nn) + (N - l)(n)] ,(35) 



00' 0" 



^ C/^(k)((2n Q -l)4c k/3 ) = (N -1)(ZU) 2 (1-2k) J ^/(£)Im{A(£)G 3 (£)}. 



(36) 



k0 

Introducing the notations 



J ^/(OIm{A(OG fl (0} = AG fl (37) 

we, for example, obtain the following formuale 

= [(N-l)(N-2)(nn) + {N-l)(n)-(N-l) 2 (n) 2 ]U 2 , (38) 

= (N - 1)(N -2)(N- 3)(nnn)U 3 + (N — 1)(N - 2)(nn)(3U 3 + e f U 2 - 2(N - l)(n)U 3 ), 

+ (N - l) 3 (n) 3 U 3 -(N- l) 2 (n) 2 (2U 3 + e f U 2 ) + (N — l)(n)(U 3 + e f U 2 ) + (N — 1){ZU) 2 {1 - 2k)AG 3 (39) 

VII. APPENDIX B 

In the crystal field case we assume that TV-fold degenerate impurity level e/ is split by a crystal field onto G sublevels 
e/i, ...£/ Q , ...e/G- We assume that for each sublevel there is still some partial degeneracy d a so that X^ a =i ^« = N. 
In limiting case of SU(N) degeneracy, G = l,d\ = N, and in non-degenerate case, G = N,d\ = d a = da = 1- We 
need to discuss how a number of electrons n can be accommodated over different sublevels e f a . Introducing numbers 
of electrons on each sublevel, n a , we obtain J2a=i n a = n - Note the restrictions: < n < N, and < n a < d a . In 
SU(N) case, G = 1, ni = n, and in non-degenerate case, G = N , n a is either or 1. Total energy for the shell with 
n electrons depends on particular configuration {n a } 

G 1 

E ni ...n G = E e f aUa + ^U(T, a n a )[(i: a n a ) - 1]. (40) 

Many-body wave function is also characterized by a set of numbers {n Q }, i.e. \n\...nQ). Energy E n ^...n G remains 
degenerate, which can be calculated as product of how many combinations exists to accommodate electrons in each 
sublevel, i.e. G^J x ...C^° x ...G^= . Let us further introduce probabilities ip ni ... nG to find a shell in a given state with 
energy £ I 7H ...„ G . Sum of all probabilities should be equal to 1, i.e. 

di d a d G 

J2 ■•• E - E Ct--Ct---Ct<...n G = !• ( 41 ) 

m— n a — riQ—O 
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There are two Green functions in Gutzwiller method: impurity Green function Gf(ioj) and quasiparticle Green 
function G g (iu>) = b Gf(iu>)b , where matrix coefficients b represent generalized mass renormalizations parameters. 
All matrices assumed to be diagonal and having diagonal elements numerated as follows: Gi(iui), ...G a (iu/), ...Gg(«w). 
Each element in the Green function is represented as follows 

VjJ — \ a — b^A a (tuj) 

Gfa{iio) — b 2 a G ga {iuj). (43) 
and determines a mean number of electrons in each sublevel 

n a = d a TY,G ga {iuy^ + . (44) 

The total mean number of electrons is thus: n — Yl a =i Hybridization function A(iu>) is the matrix assumed to 
be diagonal and having diagonal elements numerated as follows: Ai(icj), ...A a (iu>), ...Ac(iu))- Mass renormalizations 
Za — b 2 a arc determined in each sublevel. 
Diagonal elements for the self-energy are 

T, a (iu) =iui- e fa - A a (iu;) - G^iiu;) = iu(l - - e fa - (45) 



Here: 



d\ d a d G 

b a = R a L a ... ... C^...C^_\...C^'lp ni ... na ... nG 'lpn 1 ...n ct -l...n G , (46) 

ni— n a — 1 n G =0 

/ di d a d G \ - 1 / 2 

L <* = i- E - E - E ct-<:-l--c%<... na ... nG , (47) 

\ ni—O n a —l n G =0 / 

I di d a -l d G \-l/2 

R « = i- E - E - E ^■■■Cnr 1 -c^<...„ Q ...„ G . (48) 

\ ni— n Q — nc-0 / 

The generalization of the non-linear equations (|14|l has the form 
= [E ni ... nG +A-(T,®\ a n a )]-ip ni ... nc! + 

G 

{iw)G ga {iu>)]b a [Ra.La4>ni...n a -l...nG + b a 

G 

ni ■ ...tig ] ■ (49) 
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